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ABSTRACT 

We have measured the stellar velocity dispersions (a*) and estimated the cen- 
tral black hole (BH) masses for over 900 broad-line active galactic nuclei (AGNs) 
observed with the Sloan Digital Sky Survey. The sample includes objects which 
have redshifts up to z = 0.452, high quality spectra, and host galaxy spectra 
dominated by an early-type (bulge) component. The AGN and host galaxy spec- 
tral components were decomposed using an eigenspectrum technique. The BH 
masses (Mbh) were estimated from the AGN broad-line widths, and the veloc- 
ity dispersions were measured from the stellar absorption spectra of the host 
galaxies. The range of black hole masses covered by the sample is approximately 
10 6 < Mbh < 1O 9 M . The host galaxy luminosity- velocity dispersion relation- 
ship follows the well-known Faber- Jackson relation for early-type galaxies, with a 
power-law slope 4.33 ± 0.21. The estimated BH masses are correlated with both 
the host luminosities (L H ) and the stellar velocity dispersions (a*), similar to 
the relationships found for low-redshift, bulge-dominated galaxies. The intrinsic 
scatter in the correlations are large (~ 0.4 dex), but the very large sample size 
allows tight constraints to be placed on the mean relationships: Mbh oc £^ 73±a05 
and Mbh oc o-3.34±o.24 am p^ U( j e Q f the Mbh — o"* relation depends on the 

estimated Eddington ratio, such that objects with larger Eddington ratios have 
smaller black hole masses than expected at a given velocity dispersion. While this 
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dependence is probably caused at least in part by sample selection effects, it can 
account for the intrinsic scatter in the Mbh — c* relation, and may tie together 
the accretion rate with physical properties of the host bulge component. We find 
no significant evolution in the Mbh — cr * relation with redshift, up to z ~ 0.4, 
after controlling for possible dependencies on other variables. Interested readers 
can contact the authors to obtain the eigenspectrum decomposition coefficients 
of our objects. 

Subject headings: galaxies: active - - galaxies: bulges - - galaxies: nuclei - 
quasars: general 



Introduction 



It is now widely accepted that all galaxies with a massive bulge component contain a cen- 
tral massive black hole (BH) . Application of stellar dynamical and gas dynamical techniques 
for measuring masses of central BHs in nor mal galaxies has led to the identification of corre 



lation s between the BHs and host galaxies (IKormendy fc Richstone I Il995l ; iMagorrian et al. 



19981 ). The tightest of these relationships are the correlations of BH mass (Mbh) with the 



galaxy bulge luminosity and the bulge stellar velo city dispersion (cr*) (IFerrarese fc Merritt 
2000l : iGebhardt et alJEoOOat iTremaine et~al1l2002h . 



Massive BHs have also been postulated in quasars and active galaxies (jLynden-bell 



19691 : iRees I Il984l ). The question of whether AGNs follow a similar BH-bulge relation as 



normal galaxies is a very interesting one, since it may elucidate the connection between the 
host galaxy and the active nucleus. Comparisons between the bulge and BH properties in 
quasars and bright Seyfert galaxies is observationally very difficult, however, since the stellar 
component near the BH is easily lost in the glare of the active nucleus. 

Reverberation mapping (e.g. Blandford & McKee 1982; Peterson 1993; Netzer & Peter- 
son 1997; Peterson et al. 2004 and references therein) of the broad emission lines in AGNs 
has been used to determine the size of the broad line region (BLR), -Rblr = ct, where r is 
the time lag between continuum and emission line variations. Assuming that the widths of 
permitted emission lines (e.g., H/3) are due to virialized gas motions in the BH potential, the 
BH mass can be estimated from the velocity dispersion AV and the radius at the location 
of the gas -Rblr- The central BH mass is then given by 



Mi 



BH 



/ 



RblrAV 2 
G 



where / is a factor of order unity that depends on the structure and geometry of the BLR. 
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Multiple line obse rvations show an anticorrelation between emission-line lags and line 
width s, i.e., r oc AV~ 2 JPeterson fc Wandel ll99gL I2OO0I : lonken fc Peterson kooilKollatschnv 
20031 ) . which strongly supports the interpretation that the BLR dynamics are dominated by 
the BH gravitational potential. Thus, reverberation mapping provides an important tech- 
nique for estimating the ce ntral BH masses of AG Ns. The most complete reverberation-based 
BH masses can be found in lPeterson et al.l (120041 ) . Another important result of reverberation 
mapping is the discovery of a simple p ower law relationship between the continuu m luminos- 
ity and the size of BLR, R BLR oc L a fcaspi et aljboool I2OO5I : feentz et alJbood ). The value 
of a using optical wavelengths and the broad Balmer emission lines, and after correcting for 
the contribution of host galax y light, is found to be ~ 0.5, which is consistent with simple 
photoionization expectations fjBentz et al.ll2006l ). The radius- luminosity relationship pro- 
vides a secondary method of estimating the BH masses, suitable for use with single-epoch 
spectra. The single-epoch spectra technique has been applied to large samples of broad 



line AGNs to estimate BH masses (IVestergaard 



Warner et al. 112003 : IWu et al. 112004 : IVestergaard 



2002 



Mclure fc Jarvis 



2004 



2002 



Netzer 



Vestergaard fc Peterson II2006I ). 



2003 



The difficulty of measuring the host bulge properties remains, however. Image decom- 
position techniques can separate the host from the central AGNs in some cases (e.g. Mclure 
et al. 2000; Peng et al. 2002, Kulbroadt et al. 2004, Bentz et al. 2006), but they do 
not provide the necessary spectroscopic information, such as the stellar absorption lines 
that are necessary for measuring host galaxy stellar velocity dispersions (er*). Host galaxy 



absorption ] 


ines of the Call triplet at a wavelength of A ~ 8550A (e.g. 


Ferrarese et al. 


2001; 


Barth et al. 


2002: 


Nelson et al. 


2004; 


Onken et al. 


2004; 


Greene & Ho 


2006b|). However, the 



Call triplet is difficult to observe for large samples of AGNs, for the following three rea- 
sons. First, the lines are often redshifted out of the observable spectral range. Second, the 
lines are usually contaminated by night sky emission and absorption features. Finally, the 
host galaxy is usually very faint as compared to the active nucleus, causing the absorption 
lines to be very weak in contrast to the nuclear continuum. When the Call triplet is not 
available, the [Olll]A5007 narrow emis sion line is sometimes used as a surrogate for esti- 



2004 



mating the stellar velocity d ispersion (INelson I l2000l : IShields et al.l 120031 ; iGrupe fc Mathur 



Salviander et al.l 120061 ) . However, while there is a correlation between the width of 
[O ill] and cr* measured from absorption lines, the scatter in the relationship is substantial 
kelson & WhittlJl996l : bnken et al-kooi botte et alihood : IWoo et alJlioOq ) . So while the 
width of [O ill] may be used as a surrogate for a * in a statistical sense, the BH mass estimated 



from [O hi] can be u ncertain by a factor of 5 (INelson I l2000l : iBoroson 1 120031 ; iBonning et al. 



20051 ; lBorosonll2005l ). Thus, the published measurements of host galaxy velocity dispersions 
for Seyfert 1 galaxies and quasars are few in number (< 100), which prevent statistically 



-4- 



significant investigations of the Mbh-c* relationship. 



Principal comp onent analysis (P CA) has been performed on both galaxy (jYip et al. 



2004af) and quas ar (jYip et al.l l2004bl ) samples from the Sloan Digital Sky Survey (SDSS; 
York et al.l |2000| ) . These studies have shown that galaxies and quasars can be classified 
based on only a few eigencoefficients in each case. The galaxy spectroscopic sample variance 
is strongly concentrated in just the first f ew eigenspectra, with over 98% of the information 
contained in the first three eigenspectra (jYip et al.ll2004al ). The quasar information is not 
as strongly concentrated as in the case of g alaxies, but the first 10 quasar eigenspectra 
account for « 92% of the sample variance (jYip et al.l l2004bl ) . At low luminosities, the 
quasar second e igenspectrum has a s trong galactic component that resembles the first galaxy 
eigenspectrum. IVanden Berk et al.l (J2006|, hereafter VB06) used separate sets of galaxy and 
quasar eigenspectra to efficiently and reliably separate the AGNs and host spectroscopic 
components. Their tests showed that the technique accurately reproduces the host galaxy 
spectrum, its contributing flux fraction, and its classification. Inspired by the success of the 
spectral decomposition, we have separated the host galaxies from the broad line AGNs in a 
large sample of spectra using the technique of VB06. Among the merits of this method is that 
a large sample of host galaxy spectra is available from which the stellar velocity dispersions 
can be measured by the stellar absorption lines. Another advantage of this method is that the 
BH masses are more accurately measured by using the host-subtracted AGNs t han by using 



comp osite spectra, especially when host galaxy contributions are not negligible (jBentz et al. 



20061 ). The purpose of this work is to study the BH-bulge relationships for active galaxies in 
a more statistically meaningful way, given that we can obtain both host and AGN properties 
for a large number of objects. 

The structure of this paper is as follows. In § 2 we describe our sample selection, and 
in § 3 we describe the data analysis. The main results of the AGN and host parameter 
correlations are included in §4. A discussion and summary follow in §5. Throughout the 
paper we assume a cosmology consistent with recent results from the WMAP experiment 
dSpergel et al1l2006h : Q m =0.3, tt A = 0.7 and H = 70 km s" 1 Mpc^ 1 . 



2. Observations and Sample definition 



The broad-line AGNs used for this study were selected from the SDSS. The SDSS is a 
project to image approximately 10 4 deg 2 of sky, in five broad photometric bands (u, g, r, i, z) 
to a depth of r ~ 2 3, and to obtain spec t ra of 10 6 galaxies and 10 5 quasars s elected from 



the imaging suryev (jFukugita et al.l Il996l ; iHogg et al.l 1200 It llvezic et al.l 120021 ; ISmith et al. 



20021 ; iPier et al.l 120031 ). Imaging observations are made with a dedicated 2.5 m telescope 
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( IGunn et al.ll2006l ). using a large mosaic CCD camera (IGunn et al.lll998l ) in a drift-scanning 
mode. The SDSS spectra are obtained using a pair of multi-objec t spectrographs tha t 
simultaneously accept 320 optical fibers each (see lStoughton et al.ll2002l ; lBlanton et al.ll2003l ). 
Each fiber subtends a 3" diameter on the sky (see § 13.61 for a discussion of fiber size effects). 
The wavelength range of each spectrum is approximately 3800-9200A, with a resolution of 
A/ A A ~ 1850. The total spectral integration time is at least 45 minutes. 

We define broad-line AGNs here to mean any extragalactic object with at least one 
spectroscopic emission line with a FWHM of at least 1000 km s -1 , regardless of the object 
luminosity or morphology. Broad-line AGN candidates are select ed from SDSS color space, 
or as unresolved matc hes to sources in the F IRST radio catalog (IBecker. White. &: Helfand 
19951 ). as described by iRichards et al.l (120021 ) . Quasars are also often identified because the 
objects were targeted for spe ctroscopy by non-quasar selection algorithms , such as optical 
matches to ROSAT sources (|Stoughton et al.l 120021 ; Anderson et al.l 120031). various c lasses 
of stars (IStoughton et al.l 120021) . so-called serendipity objects (IStoughton et al.l 120021 ) . and 



galaxies ( Strauss et al. 20021 ). The completeness of the SDSS qu asar selection algorithm i s 



close to 95% up to the i band limiting survey magnitude of 19.1 (jVanden Berk et al.ll2005l ). 



The set of AGNs studied here were selected from the list given by VB06. In that study, 
th e AGNs were drawn fr om the catalog of SDSS Data Release Three quasars described 
by ISchneider et al.l (120051 ). and from an extension to that catalog that was constructed by 
including objects with ab s olute magnitudes fainter than Mj = —22, which is the limit im- 
posed by ISchneider et al.l (120051 ) . From that sample, VB06 selected 4666 AGNs for which 
a host galaxy spectral component could be reliably decomposed, using their eigenspectrum 
technique. Those AGNs have redshifts of z < 0.752, and host galaxy fractional flux con- 
tributions, Fjj, of greater than 10% between 4160 and 4210A. The selected spectra were 
cor rected for for e groun d Galactic extinction using the Milky W ay extinction cu r ve de scribed 
by iFitzpatrick I (119991 ). and the reddening maps provided by ISchlegel et al.l (119981 ). The 
decomposition technique uses separate sets of five galaxy eigenspectra and ten quasar eigen- 
spectra to efficiently and reliably separate the AGN and host spectroscopic components. The 
technique accurately reproduces the host galaxy spectrum, its contributing fraction, and its 
classification. The details of the AGN and host galaxy decomposition are given by VB06, 
along with the decomposition parameters and various derived parameters for each object. 
For clarity, throughout this paper we define the original spectra before decomposition as the 
active galaxy spectra, the host galaxy components reconstructed with galaxy eigenspectra 
as the reconstructed host galaxy spectra, and the AGN spectra reconstructed with quasar 
eigenspectra as the reconstructed AGN spectra. 



For the current study we made additional restrictions on the initial set of 4666 AGNs. 
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Each of the decomposed components (AGN and host galaxy) must have a sufficiently high 
signal-to-noise ratio (S/N) to reliably estimate the black hole masses and stellar velocity 
dispersions, respectively Because the M B h — relation applies only to the bulge component 
of the host galaxies, we also required that the reconstructed host galaxy components be 
dominated by bulge or "early-type" spectra. Quantitatively, the following criteria were 
applied: 1) The host flux fraction Fjj, must be between 0.2 and 0.8, which guarantees that the 
AGN and host galaxy are both significant contributors to the active galaxy spectrum. 2) The 
host galaxy cla ssification angle d > H , defined using the values of the first two eigencoefficients 



(see VB06 and lYip et al.l (l2004al )). must be in the range 0° < 4>h < 15°, which selects bulge- 
dominated galaxy spectra. The spectrum classification is based on the flux within SDSS 
fiber apertures (3"). For a galaxy dominated by a bulge component in the SDSS spectrum, 
the morphology of the entire galaxy image, extending beyond the 3" aperture, may or may 
not be classified as early type. Thus, this criterion does not guarantee that the host galaxy 
would be classified morphologically as a bulge-dominated galaxy. However, this spectroscopic 
study focuses on the bulge component of the galaxies, regardless of the extended morphology. 
Within the 3" fiber aperture, the selected sample galaxies are all bulge dominated. 3) The 
active galaxy spectrum must have a mean S/N per pixel > 15, averaged over the SDSS % band 
sensitivity function. We also required the redshifted position of the H/3 line to be covered 
by the spectra, but that constraint was already satisfied by the z < 0.752 limit in the VB06 
sample. In practice the S/N requirement limited the maximum redshift to be far below 0.75. 
These criteria resulted in a sample of 960 spectra, with redshifts from 0.013 to 0.489 (the 
median redshift is 0.153) and S/N from 15.0 to 65.7 (the median S/N is 21.6). Figured] 
shows the relation between the i band S/N and the apparent magnitude mj. Figured shows 
the S/N level as a function of redshift. Figure[3] shows the AGN component luminosity L 5W0 
(the monochromatic continuum luminosity at a rest wavelength of 5100A (see § 13. 2p ). as a 
function of redshift for the sample. The correlation of luminosity with redshift is due mainly 
to the apparent magnitude limit imposed by the SDSS quasar selection algorithm. 

It is possible that some low redshift, extraordinarily luminous objects are missed because 
of the SDSS upper flux limits (i = 15). However, these should be relatively rare objects since 
the SDSS completeness to broad-line AGN is extremely high, ~ 95% (Vanden Berk et al. 
2005). Moreover, since the relative brightness of the AGN and host is one of our selection 
criteria (0.2 < Fh < 0.8), these most luminous objects would not be selected for our sample 
since the flux is dominated by nuclear light. At the other extreme, the sample certainly 
omits the least luminous nuclear sources, both because of the SDSS faint flux limits, and 
because host galaxy light would dominate the flux. These objects would presumably, though 
not certainly, have properties that lie somewhere between those of the relatively brighter 
objects in our sample, and the galaxies with quiescent supermassive black holes examined 



- 7- 



in other studies (e.g. iTremaine et al.ll2002l ). 



3. Data Analysis 



3.1. Estimation of Black Hole Mass 



The BH masses were estimated by equation (JTJ), which requires the widths of the broad 
emission lines AV^ or FWHM, the distance of the broad emission line gas from the central 
gravitational source -Rblr, and the scale factor /. The broad emission line widths can 
be measured directly from the spectra. The gas radius cannot be measured directly from 
single-epoch spectra, so we use the result based on reverberation mapping studies that 
the continuum luminosity and radius are strongly correlated. The radius-luminosity (R- 
L) relationship has been calibrated using the gas distances from samples of objects for which 
reverberation mapping techniques could be used. The power-law index for the correlation 
between gas radius and the optical luminosity at 5100A, L5100, was reported to be about 



1 

0.67-0.70 by lKaspi et al.l (120001 . 120051 ). However, the R-L relationship can be affected by host 



galaxy contami nation; using new H ubble Space Telescope (HST) imaging of reverberation- 
mapped AGNs, iBentz et al.l ( 120061 ) corrected the host galaxy contamination of the optical 
continuum luminosity and reported a flatter power-law index of 0.52. Since our study deals 
with the nuclear lumin osities, separ a ted fr om the host stellar contribution, we adopted the 
R-L relationship from IBentz et al.l (120061 ) . The value of the scale factor f is un k nown . 
However, using the reanalyzed reverberation mapping BH masses of iPeterson et al.l (120041 ) 
and assu ming that the active g alaxies obey t he sa me M B h — relationship as quiescent 
galaxies (ITremaine et al.l |2002| ) , lOnken et al.l (120041 ) derived a value for / by requiring the 
two techniques yield statistically consistent BH m asses. In this paper we use th e most 
recently determined mass scaling relationship from IVestergaard fc Peterson I (120061 ). which 
wa s calibrated b y comb ining the R-L relationship by Bentz et al. ( 20061 ) and the scale factor 
by lOnken et al.l (120041 ) . The mass scaling relationship for single epoch spectra, using the 
FWHM of the broad H (3 emission line, and the monoc hromatic continuum luminosity at 
5100A, L5100, is given as ( Vestergaard fc Peterson 1120061 ) 



log M BH = log 



FWHMh^ 2 / AL 5100 
10 3 kms"V Vl0 44 ergs 



0.5 



+ (6.91 ±0.02) 



(2) 



where Mbh is the BH mass in solar units. The zero point error in equation ([2]) is a formal 
statistical error. A more representative uncertainty in the zero point for an individual object 
is the intrinsic scatter about the relationship of ±0.43 dex. The measurements of the H/5 
FWHM and 

£5100 f° r our sample of AGNs are described in the following subsections. 
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3.2. Emission Line Width Measurement 

The 960 spectra were first decomposed into host galaxy and AGN spectra using the 
eigenspectrum method described by VB06. Examples of the spectroscopic decomposition of 
the objects with different S/N are shown in Figs. HI [5], and [61 One can also refer to Fig. 9 
in VB06 for more decomposition examples. In each of the examples, the original spectrum 
is shown along with the reconstructed AGN and host galaxy spectra and the fit residuals. 
The reconstructed spectrum is a close match to the original active galaxy spectrum, except 
for some larger residuals in the regions of narrow and broad emission lines. Accurately 
reconstructing the emission lines can require a large number of eigenspectra, which may 
introduce spurious features into the reconstructed continuum (VB06). 

To measure the AGN emission lines, we do not use the reconstructed AGN spectrum, 
but instead use the original active galaxy spectrum with the reconstructed host spectrum 
subtracted. Because we do not use the reconstructed AGN spectrum for our analysis of 
the emission lines, the imperfect modeling of the AGN broad lines as seen in the residuals 
in Fig.H] does not affect the FWHM measurements of the broad lines. In subtracting the 
host spectrum, we are very careful in dealing with the narrow lines. After decomposing an 
active galaxy spectrum, the reconstructed host galaxy spectrum is interpolated across the 
narrow line regions, effectively removing the narrow lines. The interpolated reconstructed 
host spectrum is then subtracted from the original active galaxy spectrum. The resulting 
AGN spectrum thus contains both the host and AGN components of the narrow line regions. 
This method is better at preserving the details of the AGN broad emission line profile, and 
it retains the noise of the original spectrum, which is necessary to assess the quality of 
the spectral measurements. We define the host-subtracted active galaxy spectrum as the 
decomposed AGN spectrum. The decomposed AGN spectrum should be differentiated from 
the reconstructed AGN spectrum. As mentioned above, the decomposed AGN spectrum 
has had the host spectrum subtracted, and it retains all of the necessary original spectrum 
information, such as emission lines and noise. However, the reconstructed AGN spectrum 
is reconstructed by the eigenspectra and the eigencoefficients — this spectrum is not used 
for any spectroscopic measurements. Figs.d [H and [9] show examples of decomposed AGN 
spectra, which clearly show that the emission lines are preserved, even in the narrow line 
regions. The monochromatic luminosities at 5100A of the sample objects were measured 
from the rest frame decomposed AGN spectra, by averaging over a 10.57A (10-pixel) wide 
region centered on 5 100 A. 



The H (3 profile is not always sufficiently strong in our spectra to accurately determine 
its FWHM. We a l so me asured the FWHM of Ha for all objects possible because, as shown 
by I Greene fc Hoi (120051 ) , one can use the H a line width to approximate the width of the 
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H (3 line when the latter is unavailable. This issue and the results are discussed further in 
§ 13.31 The widths of the emission lines were determined by fits to the line profiles in the 
regions near each of the lines. To fit the H a and H f3 lines, a featureless power-law continuum 
and Fe II template were first fit to the decomposed AGN spectrum. The H a and H (3 lines 
were fitted separately. We adopted a local power law to fit the underlying continuum; the 
ostensibly line-free continuum end points were selected as 4270A and 5600A for H j3 lines, 
and as 6270A and 6980A for Ha lines. The average fluxes of the end points, averaged 
over the 21 pixel wid e regions centered on th e end points, were used to fit the power law 
continuum. Following iBoroson fe Green! (119921 ) . an effective Fell template can be generated 
by simply broadening and scaling the Fe II spectrum derived from observations of the narrow 
line Seyfert 1 galaxy I Zw 1 (kindly provided by T. Boroson). Al though a newer template 
is available (jVeron-Cetty et al.l 120041 ). we still use the template of IBoroson fe Green! (1992J) 
because one of our aims is to com pare our results to those of previous studies, all of which 
used the IBoroson fe Green! (jl992j) template. The Fe II template strength i s free to vary, but 
it is b roadened to be co nsistent with the FWHM of H/3 for each object (IBoroson fe Green 
19921 ; IShang et al.l I2005T ) . In order to avoid contamination from strong spec tral lines, we 
inclu de only the following regions in the Fell fit: 4450-4750A and 5100-5350A (IGreene fe Hoi 
20051 ). The best Fell template is found by the minimization of the \ 2 values in the fit region. 
The Fell template is very weak and may be negligible in the local region of the Ha line. 
Examples of continua and Fe II template fits to decomposed AGN spectra in our sample are 
shown in Figs. El and [91 

After subtracting the continuum and iron template, the remaining broad and narrow 
emission features near H (3 and H a were fit with a series of Gaussian profiles. The broad 
component of the H/3 line was fit with two Gaussians, w hich are usually s ufficient to ac- 
count for the width and possible asymmetry of the line (IShang et al.ll2005l ). Each of the 
[O m] AA4959, 5007A lines and the narrow component of the H (3 line were fit by single Gaus- 
sians. The broad component of the Ha line was also fit by two Gaussians, and a single 
Gaussian was used for each of the [Nll]AA6548, 6583A lines and the narrow component of 
the Ha line. The relative positions of the narrow Ha and [N n] lines were constrained by their 
laboratory values, and the relative ratio of the two [Nil] components was fixed to 2.96. For 
the H (3 line fits, however, we did not set as many strict constraints as we did for H a. Instead, 
we allow a few A shift for each line position, and let the peaks of the [O III] doublet lines 
vary independently. The H (3 constraints are looser than those for H a because H (3 is weaker 
than H a. For example, if many restrictions are placed on the fit to the [O ill] doublets, the 
fit to the broad H (3 line can become unreliable in some cases. Although the peaks and fluxes 
of the [O III] doublet lines were allowed to vary independently, the doublet flux ratios in our 
sample (Fig. lTOl) are close to the theoretical value of 3:1. All of the Gaussian components 



-10- 



within either the H/? or Ha set were fit simultaneously, with the narrow-line dispersions all 
constrained to be smaller than 5A (308.6 km s^ 1 ). The best-fit parameters were determined 
through a minimization of the \ 2 statistic, using the non-linear Levenberg-Marquarrdt min- 
imization algorithm as implemented by mpfit in Using the fitted model profiles and 
excluding the narrow components, we calculated the final FWHM of the H (3 and H a lines 
based on an analytic solution from the sum of the two fitted broad Gaussian components. 
The SDSS spectrograph wavelength dispersion (160 km s" 1 ) does not increase the broad 
emission line widths by any more than 1% even in the worst case. Example fits to the H /3 
and Ha line regions are shown in Fig.[TTl which shows the decomposed AGN spectrum, the 
fit, the narrow and broad components of the fit, and the fit residuals. 

To estimate the errors in the line width measurements, we used a method similar to that 
of Greene fe Ho J2005h . An artificial spectrum was generated for each object by combining 



the best-fit model emission lines with Gaussian random errors according to the object's error 
array. Then the artificial spectra were fit using our fitting procedure, as outlined above. 
Each spectrum was simulated by 100 realizations (a larger number of realizations did not 
significantly improve the error estimation), and the estimated error was simply taken to be 
the dispersion in the 100 measurements. The error measured in this way is the formal error, 
which accounts only for the random noise error in the FWHM measurement. Additional 
uncertainties come from the uncertainty of the fitted local continuum, and imperfect removal 
of the Fe II components. Assuming a single Gaussian profile, if the continuum is lowered by 
10% of the line peak, the estimated FWHM will increase by 7%. If we do not remove the 
Fe II components at all, the typical measured FWHM changes by about 2%. 



3.3. The Relationship Between the FWHM of Ha and H/3 

Inspection of the H/3 region of the spectra showed that some of the broad H/3 lines 
are either missing (e.g. due to problems with that part of spectrum) or are not detectable 
in the spectra. In contrast, when Ha is in the observable spectral range (z < 0.36), the 
broad component is almost always easily detectable. The Ha line can be used to estimate 
the width of the H (5 line in cases where the H j3 line cannot be measured. We used the H /3 
line to estimate Mbh when the S/N of the peak of its broad component in the decomposed 
AGN spectrum is at least 3. That criterion is satisfied in 640 of the 960 objects in the full 
sample. Among the 960 objects, the Ha line is not measurable due to bad pixels in 11 
cases and because it is redshifted outside of the observed wavelength range in 17 cases. In 
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38 cases, the FWHM of the Ha line is smaller than 25A (1143 km s -1 ), so that the narrow 
and broad components could not be unambiguously separated; these spectra were excluded 
from further H a analysis. In eight cases, the profiles of the H a and H (3 lines possess double 
emission profiles and extremely broad wings, making the meaning of FWHM ambiguous in 
these cases; these spectra were also excluded from further analysis. In total, there are 886 
spectra for which the H a line width was reliably measured. 

Among the 640 good H/3 and 886 good Ha spectra, there are 597 spectra for which 
the widths of both lines can be measured. The FWHM of H a and H (3 for these spectra 
are plotted in Fig. [12] it is clear that there is a relatively tight relationship between the two 
line widths, FWHM Hq and FWHM H /3- When fitting the relationship, neither FWHM Ha nor 
FWHMh/j can be regarded as the independent variable, so we calculate the best fit line to the 
logarithmic quantities using the "BCES" (Bivariate Correlat ed Error and intrinsic Scatter) 
bisector estimator described by lAkritas fc Bershadyl (119961 ). which takes uncertainties in 
both parameters into account. The best fit relationship using the BCES bisector estimator 
is 

FWHM H « , s /FWHM Ha \ -" ±0 - 02 i 



This result is similar to that found by I Greene fc Hoi (120051 ). and the scatter of equation ([3]) 
is about 0.06dex (~ 470 km s _1 ). In Fig.[12l the solid line shows the best fit, and the 
dash ed line denot e s FW HMh„= FWHMh/3- The filled dots are the mean spectrum FWHM 
from iKaspi et aJj (120001 ) for comparison; those points are also described well by the best fit 
relation. 



3.4. Uncertainty in BH Mass Estimates 

The uncertainties of the FWHM measurements for H a and H (3 and of the continuum 
luminosity include not only the errors from the fitting routine, but also the errors that arise 
in the decomposition of the AGN spectrum from the original active galaxy spectrum. We 
estimated the total uncertainty by using simulations in the following way. For a typical object 
with the median S/N of the sample, we simulated 100 spectra by adding random noise to 
the original active galaxy spectrum according to its error array. Each of the 100 simulated 
spectra was decomposed, and the continuum luminosity and FWHM of the emission lines 
were fit as described above. 

The uncertainty in each parameter was calculated by measuring the rms dispersion of 
the 100 measurements. The typical rms dispersions for the FWHM measurement of Ha 
and H (3 are 4% and 7% respectively; the luminosity at 5100A has a typical dispersion of 
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4%. The final uncertainty for the BH mass by equation (j2J) is 7% and 14% respectively, 
using the FWHM of Ha and H/3. A representative value of the uncertainty in the BH 
mass for an individual object is the zero point uncertainty of the mass scaling relationship, 
e.g. equation (|2l). The standar d deviation of the zero point (intrinsic scatter) is ±0.43dex 
( jVestergaard fc Peterson II2006I ). 



FigurefTBl shows the distribution of the AGN monochromatic continuum luminosity, 
£5100 • The solid line is the distribution of L 5100 after subtracting the host galaxy component. 
The dotted line shows the estimated L 5W0 values when the host galaxy component is not 
subtracted; these values were estimated directly from the extinction corrected active galaxy 
spectrum at rest-frame 5100A, averaged over a 10-pixel wide region. It is clear that the host 
galaxy component is non-negligible in our sample; the difference between the luminosity with 
and without host galaxy subtraction is close to 0.5 dex. This implies that the host galaxy 
components can significantly affect the accuracy of the BH mass estimates at the luminosities 
present in our sample. The host galaxy contri bution will also affect the luminosity-radius 
relationship, as discussed by lBentz et al.l (120061 ). 



3.5. BH Mass Distribution 



The black hole masses were calculated according to equation ([5]). The FWHM of the 
H p line was used directly in the 640 cases for which the line could be measured. For the 286 
cases in which the H a line was measurable when H (3 was not, the H a FWHM was used to 
estimate the FWHM of the H (3 line according to equation (j3J). The distribution of black hole 
masses is shown by the dotted histogram in Fig.O Most of the black hole masses are in the 
range 10 6 < Mbh < 1O 9 M . From the black hole masses and the continuum luminosities, we 
can estimate the Eddington ratio, Lboi/-^Edd — the ratio of the AGN bolometric luminosity 
to the Eddington lu minosity for the bla ck hole — using the approximate relation Lboi ~ 
10AL A (5100A) (e.g.. IWandel et al.lll999l ). The estimated Edd ington ratio has a st andard 



deviation of 25% due to variations in the bolometric correction (IRichards et al.ll2006l ). Since 
-^Edd is proportional to Mbh, the Eddington ratio would have a standard deviation of ~ 28%, 
assuming the typical error in our BH mass estimates (§ 13. 4p . The solid histogram in Fig.[T4l 
shows the distribution of -kboi/^Edd; the Eddington ratios also fall within a limited range, 
mostly at small values. The small Eddington ratios and limited range of BH masses are 
due to the sample selection criteria, which require that the stellar absorption features are 
detectable at a level that allows the host galaxy velocity dispersions to be measured. 

Fig. [15] shows the relationship between BH mass and the estimated Eddington ratio 
Lboi/^Edd- I n our data sample the BH mass is inversely correlated with the Eddington ratio. 
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This is not unexpected if the range of bolometric luminosities does not change proportion- 
ally with Mbh- Both Mbh and Lb i/-^Edd are parameters calculated from AGN continuum 
luminosity Lagn ~ AL A (5100 A), and emission line FWHM. Therefore, Fig.[T5] is simply a 
remapping of the sample L^gn and FWHM values. The correlation is due in part to sample 
selection effects in the SDSS and our sample. Since our sample is selected only from SDSS 
AGN with detectable stellar features (0.2 < Fh < 0.8), the luminosity cannot be above 
the SDSS bright limit for spectroscopy (i = 15). Our sample is neither so bright that the 
flux is dominated by nuclear light, nor so faint that the flux is dominated by host galaxy 
light. Thus, objects with higher BH masses tend to have a lower range of Eddington ratios, 
and vice versa. The dispersion in the relationship is largely due to the range of measured 
emission line widths that are used to calculate Mbh- 



3.6. Host Galaxy Velocity Dispersion Measurement 



Two accurate and objective methods, the Fourier-fitting method and th e direct-fitting 



method, have been developed for measuring the stellar velocity dispersion (jSargent et al. 



1977 



2002 



Tonrv &: Dayis 



19791 : iFranx et al.lll989l : iBender lfl99~ol ; |Rix fc Whitelll992l ; iBarth et al. 



Greene fc Holl2006bl ) . They are all based on a comparison between broadened template 



spectra and the spectrum of the galaxy whose velocity dispersion is to be determined. Fourier 
space is the natural choice to estimate the velocity dispersions because a galaxy's spectrum 
is a mix of stellar spectra convolved with the distribution of velocities within the galaxy. 
However, there are several advantages to treating the problem entirely in pixel space, rather 
than Fourier space. Specifically, the effects of noise are much more easily incorporated in 
the pixel-space-bas e d direct-fitting method . We have used the "direct-fitting" method (e.g. 
Rix fc White! Il992l ; iBernardi et al.l l2003al ). in which the spectrum is directly fit in pixel 
space. We use the public IDL program vdispfit, written by David Schlegel, to find the 
velocity dispersions for the host galaxies. The templates consist of the first four eigenspectra 
from a principal comp onent analysis of the echelle stellar spectra in the Elodie database 
(IMoultaka et al.l 12004 ) . The best-fitting dispersion value was determined by minimizing \ 2 
for the fit. 

The reconstructed host galaxy spectrum was not used to measure the velocity disper- 
sion, because features such as narrow absorption line profiles may not be well reconstructed 
without invoking a much larger number of galaxy eigenspectra than we have used here. In- 
stead, we used the observed active galaxy spectrum after subtracting the reconstructed AGN 
component, which we define as the decomposed host galaxy spectrum. This method allows 
us to use the error array of the original spectrum to determine measurement uncertainties. 
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In the fitting routine, all of the narrow emission lines were blocked, because there is very 
little information about the absorption spectrum in those regions, and the AGN fit to the 
narrow emission lines may be poor. We also blocked the broad emission line regions of H a 
and H f3 to avoid any artifacts introduced by imperfectly reconstructed AGN components. 

An example fit of a host galaxy (from among the same objects used as examples in 
Fig. ED, resulting in a measurement of the velocity dispersion, is shown in Fig.[16j In the 
figure, the two large boxes show the masked broad Ha and H/3 line regions; the left and 
the right small boxes show the masked narrow emission line regions (H7 and [Oi]), and the 
middle small box shows the masked bad pixel region. 

The SDSS spectra were obtained by using a fixed 3" fiber, which covers different fractions 
of the projected areas of different objects, due both to different intrinsic galaxy sizes and dif- 
ferent redshifts. Thus, the velocity dispersions should not be compared without accounting 
for the finite fiber diameter, which is discussed here. The dispersion is scaled to a stan- 



dard relative circular aperture, defined t o be one-eighth o: 



correction formula was applied following Ijprgensen et al.l (119951 ) and IWegner et al.l (119991 ) 



the effect ive galaxy rad i us. T he 



°cst 



ffibe: 

r /l 



0.0 i 



(4) 



where r^r = 1"5 an d r is the effective radius of the galaxy measured in arcseconds. 
Bernardi et al.l (j2003al ) studied early type galaxies up to a redshift of z = 0.3 and showed 



that most gala x ies hay e r > l'/5. Our sample has almost the same redshift limit as that of 
Bernardi et al.l (j2003al ) and the correction depends only weakly on r , so the correction is 
not large. The effective radii of the galaxies in our sample have not yet been measured, so we 
assume that all of the early type galaxies (or the galaxy bulges) in our sample have typical 
effective radii, which corresponds to an angular size of Tq = l'/5 at a redshift of z = 0.3. At 
redshifts other than z = 0.3, the angular sizes were scaled with distance. The correction has 
a maximum value of 10.8% and a median value of 5.4%. For a more accurate correction, the 
effective radii of the hosts must be measured. In the case of inactive galaxies, estimations 
of effective radii are publicly available parameters measured from SDSS images. For active 
galaxies, one should first remove the central nucleus to get the host galaxy effective radius, 
which requires image decomposition to separate the nucleus and host. We have not pur- 
sued this additional correction here, owing to the limited improvements possible with SDSS 
imaging data in relation to the required workload. 
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3.7. Uncertainty in the Velocity Dispersion Measurement 



It was shown by lBernardi et al.l (l2003al ) that the direct-fitting method does not produce 
large systematic errors, but the same conclusion may not apply to our measurements given 
the presence of AGN components. Even though a measured velocity dispersion may be 
accurate for a decomposed host galaxy spectrum, it may be different from the true value if 
the decomposed host galaxy is not an accurate representation of the true host galaxy. A 
test for systematic errors in the velocity dispersion measurement was made by generating 
simulated data. We first constructed an early type galaxy template using the stellar echelle 
spectra in the Elodie database. The template early type galaxy has a high S/N and a very 
small velocity dispersion. The template was broadened with different values of the velocity 
dispersion, ranging from 60 to 400kms _1 , and rebinned to the same wavelength scale as 
the SDSS spectra (69kms _1 /pixel). To each of the broadened galaxy templates, a template 
AGN spectrum was added to form simulated active galaxy spectra. In all cases, the fractional 
contribution of the host galaxy was set to F H = 50%. 

The simulated active galaxy spectra were decomposed into host galaxy and AGN com- 
ponents, and the direct-fitting method was used to make measurements of the velocity dis- 
persions. No noise was added to the spectra at this point, so any deviations of the measured 
velocity dispersions from the true values represent systematic errors in the measurements. 
Figured?] (dotted line) shows the relative difference between the true velocity dispersions and 
the measured values, 5a* = (o-*^ measure d — o~*jrue) / &*,true- The measured velocity dispersion is 
systematically approximately 1% smaller than the true velocity dispersion. This systematic 
error steadily increases to about 4% below velocity dispersions of about 120 km s -1 . 

To test the effects of noise, we added random noise to the template active galaxy spectra. 
In our sample, the lowest spectroscopic S/N is 15, which sets the worst case limit in the 
simulations (of course, objects with a host galaxy fraction of less than the simulation value 
of F H = 50% will usually return poorer results; see below). Random noise corresponding 
to a 5/iV of 15 was added to each template 100 times. The noisy spectra were decomposed 
and the velocity dispersions were measured as usual. The dashed line in Fig.[T7J shows the 
relative error between the true velocity dispersion and the mean measurement in each set 
of 100 simulated spectra. The solid curves in Fig. [17] show the rms scatter of the velocity 
measurements from the mean, due to the random error in 100 simulations. Fig.[T7J shows 
that the random error of the velocity dispersion measurement is around 10% for the spectra 
with the lowest S/N in our sample. The mean relative errors track the noiseless relative 
errors very closely, showing again an approximately 1% systematic offset. The systematic 
relative error is larger at smaller velocity dispersions, reaching a maximum absolute value of 
about 4% at 60kms _1 , which is close to the dispersion per pixel of the SDSS spectrographs 
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(69 km s l ). Values of the velocity dispersions smaller than this limit are very uncertain. 

The vdispfit routine returns an estimate of the uncertainty in the velocity dispersion, 
based on the spectral noise, which we find to be consistent with the results of the simulations. 
Figure[18] shows the estimated uncertainty in the measurement of the velocity dispersion, 
returned from the vdispfit routine, as a function of the decomposed host galaxy S/N. As 
expected, the errors are generally smaller with increasing S/N. Some of the errors are 
greater than the worst case results from the simulations, because their host galaxy fractions 
are less than F H = 50%. FigurefTH] is the simulation of the the worst case (F H = 20% and 
S/N = 15), which shows that the maximum error can be 30% for small velocity dispersions. 
Error may arise if the templates are not good matches to the host galaxy spectra. For 
example, the templates may not be able to account for a strong post-starburst component. 
However, all of our host galaxies are chosen so as to be dominated by old stars. In addition, 
the simulations did not account for bad pixels, night sky line residuals, and other artifacts. 
The reduced \ 2 values in the velocity dispersion fittings are close to 1, which proves that 
any template mismatch is negligible in our sample. The vast majority of the objects have 
velocity dispersion errors of less than 15%. 



3.8. The AGN and Host Galaxy Data Set 

Spectra whose velocity dispersions are less than 69kms _1 were excluded from further 
analysis. After all of the selection criteria were applied, there were 617 objects with reliable 
H/? and 840 with reliable Ha line measurements, all of which also have reliable velocity 
dispersion measurements. 

Black hole masses were calculated using the B./3 FWHM in the 617 cases for which it 
was measurable, and in another 286 cases, using an estimate of the H/3 FWHM based on 
the Ha measurements and equation ([3]). The total number of objects (out of the initial 960) 
with reliable parameter measurements is 903. The electronic version of Table [T] lists all of 
the active galaxies and their measured or derived parameters, including redshift, L 5100 , Ha 
and H/? FWHM, Mbh, ct*, and host luminosity (as calculated by VB06). A partial table is 
shown in the printed version of Table [TJ 
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Results 



4.1. The fitting algorithm 

In the following sections, we will study how the BH masses, host galaxy luminosities and 
the host velocity dispersions are correlated. We assume that there is an underlying linear 
relation in the logarithmic quantities for each pair in the form y — a + f3 x. The best fit 



param eters a and (3 were found by the minimization of x 2 , defined, following iTremaine et al. 
fl2002h . as 



X 



N 

E 



a 



f3 Xi f 



e 2 - + (3 2 e 2 - 



(5) 



where X{, yi correspond to the measurements for each object, and e« is the formal uncertainty 
in the measurement. The 1 a uncertainties in a and (3 are given by the maximum range of a 
and f3 for which x 2 ~ Xmm < 1- To minimize \ 2 i we use t ne Levenberg-Marquardt algorithm 
as impl emented by mp fit in IDL, which recovers the results of the Numerical Recipes routine 
fitexy (IPress et al.lll992l ) as implemented in IDL for the case of symmetric errors. 



4.2. Host Galaxy Bulge Luminosity and Velocity Dispersion Relation 

The luminosities and velocity dispersions of early type galaxies have long been known 
to be correlated — a phenomenon known as the Faber- Jackson relation (IFaber fc Jackson 



19761 ). The host galaxy absolute magnitude in the SDSS g band, M g (calculated by VB06) is 
shown as a function of velocity dispersion in Fig.[20j The host galax y M g range of ou r sampl e 
is very similar to that of the normal elliptical galaxies studied by iBernardi et all (j2003al ). 
except ours are a bit more luminous. There is a wide intrinsic scatter in the Faber- Jackson 
relation, which we also find in our sample, so for clarity we have b inned the values of th e 
host galaxy luminosity in the figure (similar to the procedure used by IBernardi et al.ll2003bl ). 
Each point represents the average value for at least 50 objects. The bars in the figure show 
the rms scatter in the velocity dispersion in each bin (the errors on the mean values would be 
approximately v^50 times smaller). There is clearly a strong correlation between M g and <x*, 
although we find a sligh tly larger dispersion ab out the relation relative to studies of inactive 



early type galaxies (e.g. IBernardi et al.ll2003bl ). most likely due to the typically smaller S/N 



in our decomposed host spectra. The solid line shows the best power-law fit to the relation 



L oc a. 4 - 33±a21 



(6) 



which is consistent with the Faber- Jackson relation for inactive early type galaxies (IBernardi et al 
2003bl ). Table [2]lists the mean values of the host luminosity, stellar velocity dispersion, black 
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hole mass, and rms dispersions for the data sets in each bin. The binned values will be used 
in the analysis of the following subsections. 



4.3. Black Hole Mass and Bulge Luminosity Relation 

Figure [2T] shows the relationship between BH mass and absolute g band magnitude M g 
of the host galaxies. Crosses show the mean value in each bin, and error bars show the rms 
scatter around the mean value. Each point represents the average value for about 50 objects. 
A linear fit using the logarithmic Mbh values gives the relation 

logM BH = -0.29(±0.02)M 3 + 1.46(±0.33) . (7) 

From the relation in Eq.[?l we can also derive the expected relationship between Mbh 
and host bulge mass Mi }U i ge by assuming that only the host bulge contributes to M g . We 
adopt a g-band mass-luminosity relation M bu i ge oc L^f^ ' 03 , consistent with the V-band re- 
sults of Magorrian et al. (1998). Substituting in, we find a relation M BH oc M^f om . This 



implies that the Mbh — M-mw re 
with the result found bv IWande 



bulge 

ationship is non-linear in our sample, and it is consistent 
How e ver, i t is in contrast to the linear relationship 



del J2002h 



found by others (e.g. iMclure fc Dunlopll200ll . |2002| ). Verification of our result will require 



more precise measurements of the bulge luminosity, since our measurements rely on spectral 
decomposition, rather than image decomposition, and we have assumed that the host com- 
ponent of the SDSS g band "cmodel" magnitude contains only the bulge flux. The "cmodel" 
magnitude uses a best fit linear combination of an exponential and de Vauc. profile, which 
gives a robust estimate of the total flux from an extended object, regardless of morphology. 
The bulge flux was estimated by subtracting an appropriately scaled nuclear component 
flux, derived from the decomposed AGN spectrum, from the "cmodel" magnitude flux, as 
described by VB06. A more precise bulge luminosity could be obtained from image decom- 
position, which directly separates out the host galaxy bulge and nuclear AGN components. 
As a complement to the spectral decomposition, image decomposition is planed for a future 
paper. 



4.4. Black Hole Mass and Velocity Dispersion Relationship 



Figure [22] shows the BH mass and a* relation for all the individ ual points, and the do t 
dashed line is the best fit to the Mbh _ c* relation for inactive galaxies (ITremaine et al.ll2002l ). 
For log a* > 2.2, objects in our sample have significantly smaller average BH masses than do 
inactive galaxies at the same value of a*. The measured values of slope f3 in the M B h oc of 
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relation cover a wide range from 3.5 t o 5.0 (IFerrarese fc Merrittll2000l ; iGebhardt et al.ll2000al ). 
As studied by iTremaine et all (120021 ). the wide uncertainties in the slope are due to the use 
of different statistical algorithms, the unknown intrinsic scatter in the BH mass estimation, 
and the uncertainties in the velocity dispersions. Our large sample helps reduce some of 
these uncertainties. 

The relationship between the black hole mass and host velocity dispersion was found 
by averaging the values of the two quantities in bins of host absolute magnitude. That is, 
the values represented by the points in each M g bin in Figures [2D] and [2T] were compared 
to find the Mbh — o"* relation, which is shown in Fig. [231 The crosses show the mean value 
in each bin and error bars are the standard deviation of the mean in each bin. Each point 
represents the mean value for about 50 objects, and the error bars are those in Figuresl2"0l 
and [2T] divided by the square root of the number of objects in each bin. The solid line is the 
best-fitting linear relation (to the logarithmic quantities) 

log M BH = 3.34(±0.24) log(a*/200km s" 1 ) + 7.92(±0.02) . (8) 



The slope found here is 1.7a flatter than the one found for inactive galaxies by lTremaine et al 



( 120021 ). who found a value of 4.02 ± 0.32. Fixing the slope to that value, the intercept of 
the best fit to our data become s 7.96±0.02, which is smaller than the value of 8.13 ± 0.06 
found by lTremaine et al.l (120021 ) by 0.17 dex, a significance of 2.7<r. (In both cases a is the 
quadrature sum of the uncertainties of both measurements.) The same Mbh — cr * relation 
as equation ([8j) can also be obtained by fitting the individual data points in Figure [221 when 
an intrinsic scatter of 0.36 dex on BH masses and 0.12 dex on the velocity dispersions are 
included. This scatter is consistent with the dispersions in each binned data set (see Table [2]). 



Figure [2H shows the same relationship as Figure [23] but with values from the literature 
added in the plot. The additional data extends well beyond the dynamic range of our values, 
particularly toward smaller values of Mbh and a*. The lite rature data and uncertainties are 
all obtained from the table given b y iGreene &: Hoi (j2006af ). The o pen squares a r e the data 



Nelson et a 



Barth et al 



estimated bylGreene &: Hoi (j2006al ). the asterisks are the data fr qmlOnken et al.l (2004) and 



( 2004j) , and the open triangles are the data from iPeterson et al.l (120051 ) and 



(120041 ) for NGC4395 and POX 52 respectively. The filled circles are the binned 
data from this paper as shown in Fig. [23] In Fig.[24lthe solid line is the best fit usin g all of 



the data. The dashed line is the fit for inactive galaxies from ITremaine et al.l (120021 ). 
The best fit to the Mbh — a * relation using all of the data is 



log M BH = 3.93(±0.10) log(a*/200kms~ 1 ) + 7.92(±0.02) . 



(9) 



This slope is entirely consistent with that measured for inactive galaxies (ITremaine et al 



-20- 



20021 ). Fits were also made with a fixed slope of 4.02 to all of the data. The fixed-slope 
intercept is smaller than the inactive galaxy value of 8.13 by 0.19 dex, a significance of 3a. 

Based on our analysis, we find that the Mbh — c* relationship for active galaxies in our 
sample is marginally flatter than that of inactive galaxies, but that it is consistent with that 
of inacti ye galaxies if the liter ature data are added. The flattening is consistent with the 
study of iGreene &: Hoi (]2006al ). The intercept of the active galaxy relationship is smaller 
than that for inactive galaxies. This holds whether or no t the literature dat a are in cluded in 
the analysis, and is consistent with iNelson et al.l (120041 ) and IGreene fc Hoi (]2006al ). though 
with uncertainties ~5 times smaller than in the former reference and ~2 times smaller than 
the latter. 

We caution that the results presented in this section depended in part on the selection of 
the mass-scaling relationship used to estimate Mbh, and that not all of the literature values 
were derived using the same relationship that we use here. In particular the intercept of 
the Mbh — 0"* relationship would be even smaller had we used the commonly adopted scale 
factor value / = 3 (which is inconsistent with the value derived empirically by lOnken et al. 



( 120041 ) ). The power-law index of the radius-luminosity relationship also affects the slope 
of the Mbh — c* relationship. A higher radius- luminosity power-law index (which would 
be inappropriate for our study since we account for host galaxy flux contamination) would 
increase the slope. 



4.5. The Dependence of the AfBH — Relation on the Eddington Ratio and 

Redshift 

The distribution of data points in the Mbh — cr* plane may be affected by many factors, 
such as the details of the BLR physics and galaxy evolution. We can begin to examine the 
dependence of the distribution of points in the plane on some of these factors, at least in a 
statistical sense, thanks to the large size of our data set. 

Figurel25l shows the Mbh — cr* relationship for two different ranges of Lb i/-^Edd- The 
AGN sample was divided at the median value of the Eddington ratio, Lboi/^Edd = 0.027. 
Objects with L ho \/L Edd values below the median are shown as filled squares, and those with 
values above the median are shown as open triangles. The median values for the low and high 
-^boi/^Edd groups are 0.015 and 0.072 respectively. Each point represents the mean value for 
at least 25 objects. There is a clear separation in the Mbh — 0"* relationship between the 
high and low Lboi/^Edd samples. For a given velocity dispersion, those objects with lower 
values of L ho \/L Edd have larger black hole masses. The slopes of the relationship in the 
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two samples are not significantly different, but the intercepts differ greatly. With the slope 
fixed at 4.02 as for inactive galaxies, the intercepts are 8.15±0.02 and 7.68±0.02 for the 
samples with low and high L ho \/L Edd respectively. The intercepts differ by about 0.47 dex 
in log(M B H) in the two groups. The intercept of the low L bol /L Edd sample is a bit higher 



than that of inactive galaxies (ITremaine et al.ll2002l ). which of course have extremely low 



-^boi/^Edd values. Generally then, in the sample studied here, BH masses are smaller relative 
to the mean relation for objects with larger Eddington ratios. 

The 0.47 dex difference of the intercepts in the two groups is statistically significant 
(> lOcx), given the precision with which they can be measured with the large samples. At 
face value, this result would tie together the accretion rate with physical properties of the host 
bulge component. However, there may be selection effects and parameter interdependencies 
that could cause the apparent difference. For example, Fig. [15] shows an anticorrelation 
between Mbh and Lboi/^Edd, but as discussed in § I3.5[ it is caused by selection effects and 
the functional interdependence of the two parameters. Therefore, to determine the reality of 
the apparent dependence of the Mbh — c* relationship on Lboi/^Edd, correlations among the 
other parameters that are involved — AGN monochromatic luminosity L51005 emission line 
FWHM, host luminosity L H , stellar velocity dispersion a„, and redshift — must be taken 
into account. To test the reality of the results, we have used partial correlation analysis to 
account for the interdependencies. 

For the partial correlation tests, we defined a new variable AMbh, 

AM bh = log(M BH ) - 3.34 log(a*/200km s^ 1 ) - 7.92 , (10) 

which is the difference, in logarithmic space, between the measured BH mass and the BH 
mass predicted by equation (|S]). Fig.[2H] shows the relationship between AM B h and Ed- 
dington ratio; there is apparently a strong anti-correlation, with a Pe arson product-moment 



correlation coefficient of —0.576. A partial correlation analysis (e.g. IWall &: Jenkins 1 120031 ) 
was performed on AMbh and Lbd/^Edd, controlling for L5100, FWHM, a*, Lh, and red- 
shift; the analysis was performed in logarithmic space, because many of the relationships 
are then more closely linear. The results are shown in Table [3J Each row of Table [3] shows 
the AMbh vs. Lboi/-^Edd partial correlation coefficient, r&M,Eratio;x, where x is one of the 
control variables, and the probability that such a coefficient value or less would occur by 
chance. Controlling for each of the other variables, there is less than a 0.01% probability 
that the AMbh anti-correlation with Eddington ratio is due to chance. Table [3] lists the cor- 
relation coefficients, r Erati0jX and r& MjX , between Eddington ratio and AM B h respectively, 
and the control variables. It is easy to show from the partial correlation results that the 
anti-correlation of AMbh with Lboi/^Edd is not due to the interdependencies of Mbh, £5100, 
and L hol /L Edd . For a fixed FWHM, both M BH and L bo i/-^Edd are positively correlated with 
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^5100) by their definitions. Thus, by equation ffTUj) . one would expect a positive correlation 
between AM B h and L bol /L Edd for a given FWHM, assuming no other dependencies. How- 
ever, the partial correlation coefficient, controlling for FWHM, is negative (Table [3D and 
still highly significant. Therefore, the anti-correlation of AM B h with L bol /L Edd is not due 
to the interdependencies of Mbh,-^5ioo, and L^ \/Ledd. 

The partial correlation tests do not determine whether the AMbh vs. Lboi/^Edd anti- 
correlation is due to sample selection effects. Figure [57] shows that there is a correlation 
between AGN monochromatic luminosity and host galaxy velocity dispersion, with a corre- 
lation coefficient of 0.569. This correlation can explain the difference between the intercepts 
in Fig. [25] For a given BH mass, objects with larger values of cr* will generally have larger 
values of L 5W0 , and therefore, larger Eddington ratios. Therefore, objects to the right of 
the regression line (having larger a* values) will generally have larger Eddington ratios, and 
vice versa. However, the correlation between L5100 and a* may be a consequence of the 
Faber- Jackson relation (§ 14.21 Fig. 120]) and a correlation between L5100 and host luminosity 
Lh, shown in Fig. [28] (cf. VB06). Assuming the Faber- Jackson relation we see is not due to 
selection effects, the reality of the Mbh — cr* dependence on the Eddington ratio depends 
upon whether or not the L 5100 — L H correlation is intrinsic, or due to selection effects. As 
discussed by VB06, at least part of the correlation is due to the restriction that both the host 
and AGN components contribute significantly to the composite spectrum; this corresponds 
to the 0.2 < Fh < 0.8 criteria described in §[2] There are clearly many AGNs in the universe 
that do not satisfy those criteria. The partial correlation coefficients of AGN monochromatic 
luminosity and host galaxy velocity dispersion are shown in Table H] Although the correla- 
tion is weaker controlling for Lh and/or redshift, it is still significant. Hence, the correlation 
of AGN monochromatic luminosity and host galaxy velocity dispersion is intrinsic, but is 
enhanced due to a sample selection effect. To summarize, the apparent M B h — depen- 
dence on L bo i/L Edd remains after accounting for several possible interdependencies, although 
it is likely that the relative luminosity criteria used to construct the spectroscopic sample 
artificially strengthens the apparent dependence. 

Selection effects similar to ones imposed on the spectroscopic sample can be avoided, 
at least pa rtially, in high-resolution imaging studies, such as those carried out wit h the 



EST (e.g. 


Flovd et al. 


2004 




Sanchez et al. 


Hamilton et al. 


2002; 


Schade et al. 


200oL In 



2003 



for example, a strong correlation was found between nuclear and host luminosities. The 
explanation is that the data are bracketed by lines of constant Eddington luminosity, ranging 
from Lboi/-^Edd ~ 0.01 to 1. The cutoff at an Eddington ratio of near unity is particularly 
sharp, possibly indicating that very few AGNs emit above their Eddington limit. The low 
Eddington ratio limit is more likely to be a selection effect, in that low-luminosity nuclei are 
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difficult to detect in high-luminosity hosts. The apparent correlation, whether intrinsic or 
a selection effect, would cause at least qualitatively the same dependence of the M B h — o"* 
relation on the Eddington ratio seen in Fig.[ 



Figure[29] shows the Mbh — cr* relationship for data in several different redshift bins. 
Our sample was divided into three subgroups by redshift: < z < 0.1, 0.1 < z < 0.2, 
and z > 0.2, shown in the figure by open triangles, filled squares, and crosses respectively. 
The data were separated into bins of host M g , as with the previous analysis; each point 
represents data from at least 25 objects. The points clearly occupy different regions along 
the best fit Mbh — c* line. That is due to the redshift-luminosity selection effect correlation 
in the sample, seen in Fig. [3} objects at higher redshifts have higher luminosities on average, 
so AGNs with larger host velocity dispersions and more massive black holes are preferentially 
selected at higher redshifts. The solid line shows the best fit to the Mbh — °"* relation from 
equation (|H]). Figure [291 shows that the intercept of the relationship doesn't change greatly 
with redshift. The points in the highest redshift bin are on average slightly above the best fit 
line, while those at lower redshifts are slightly below the line on average, suggesting a possible 
redshift evolution of the Mbh — o"* relation. With the slope fixed to the value in equation (jSJ), 
the intercepts of the three groups from low to high redshifts are 7.85±0.02, 7.88±0.02 and 
8.00±0.03 respectively. At face value, the results suggest a weak but significant dependence 
of the Mbh — relation on redshift, such that active galaxies on average have larger black 
hole masses for a given velocity dispersion at high redshift as compared to low redshift. 

As with the dependence of the Mbh - c* relation on Eddington ratio, a partial correlation 
analysis was performed to determine whether the apparent redshift dependence of A Mbh 
could be accounted for by correlations with other parameters. The correlation coefficient of 
AMbh with redshift, not accounting for other parameters, is 0.005 (see last row of Table EJ), 
which is relatively small but may or may not be significant. The results of the partial 
correlation analysis, with L 5100 , FWHM, cr*, L H , and L bol /L Edd as the control variables, 
are shown in Table [51 For each of the control variables, the remaining correlation between 
AMbh and redshift is relatively small. The significance of the remaining correlations are 
all quite high (probabilities < 0.11%), except when Lu is the control variable, in which 
case the probability of a chance correlation is about 11.8%. However, when controlling for 
all five parameters, the partial correlation coefficient is 0.059 with a random probability of 
being greater of 92.4%. Thus, there is no compelling evidence for redshift evolution in the 
Mbh — o"* relation for the sample presented here. 

The lack of evidence for redshift evolution is in c o ntrast to the results of several studies 
flTreu et alJliool I Woo et alJliooS IPeng et al.l[2006aJ lbl: Ishields et al 1 l2006h. that have found 



apparent positive correlations between the Mbh — o"* relation and redshift. 



Woo et all (j200J ) 
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found that their sample of Seyfert galaxies at z ~ 0.36 had larger b lack hole masses for 
a given a*, relative to the values expected from the local relation. IShields et al.l (120061 ) 
measured CO emission lines to infer cr* for a set of nine high-redshift AGNs, and found 
that th e host masses are u ndersized by more than an order of magnitude, given their BH 
Peng et al.l (|2006al Jbh found that z > 2 hosts in their sample a re less luminous (and 



masses. 



ther efore less massiv e ) tha n expected for a given Mbh- In contrast, IShields et al.l (120031 ) 
and Salviander et al. ( 20061 ). using O [ill] line widths as surrogates for a*, found no evidence 
for redshift evolution u p to z ~ 3. The re is also evidence from a study of high-redshift 
sub- millimeter galaxies (IBorys et al.l 120051 ) that at least some galaxies at high redshift have 
BH masses substantially smaller than expected for their host masses. The methods and 
equations used to estimate Mbh and a* were somewhat different in each study, and different 
from our own, and each sample was selected in different ways, so direct comparisons are not 
necessarily valid. Our study is the first to examine the possibility of redshift evolution for a 
single, large, homogeneously collected and analyzed sample spanning both low and relatively 
high redshifts; and our study is the first to determine whether an apparent dependence of 
the redshift evolution could be the result of correlations among other parameters. It will be 
important to extend the studies of the M B h — c* relation to higher redshifts, over a wide 
range of BH masses and velocity dispersions, using a variety of techniques. 



5. Discussion and Summary 

The very large and homogeneous data set analyzed here has allowed parameterization of 
several relations for broad-line active galaxies. Host galaxy luminosity and velocity dispersion 
are described by the Faber- Jackson relation, central black hole mass is correlated with host 
bulge luminosity and velocity dispersion, the amplitude of the M B h — relation depends 
inversely on the Eddington ratio, and there is no significant change in the Mbh — c* relation 
with redshift up to z ~ 0.4. 

The dependence of the Mbh — a * relation on both redshift and Eddington ratio have 
been matters of dispute. In our sample, we find no significant correlation with redshift 
after accounting for possible correlations among many other parameters. As discussed in 
§ 14.51 there are studies that support positive, negative, and no correlation with redshift. 
The M B h — °"* relation appears to depend on the Eddington ratio, even after accounting 
for other parameters. Apparently, the difference we have found between the high and low 
-^boi/^Edd samples is sufficient to account for the intrinsic scatter in the Mbh — <?* relation, 
usually measured to be about 0.4 dex in black hole mass for a given velocity d ispersion (e.g. 



Greene fc Holl2006al . §H]of this paper). We checked this result by dividing the lGreene &: Ho 
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( )2006al ) sample by its median Lboi/^Edd; this sample shows the same qualitative trend that 
objects with higher Eddington ratios have smaller black hole masses for a given velocity 
dispersion (although the uncertainty is too great to say whether it can account for the 
scatter). As discussed in § 14.5} the dependence of the M B h — relation o n L hn] /L FM is 
probab ly at least partly due to selection effects — which may also affect the I Greene fc Ho 
( 12006al ) sample — but a significant intrinsic component cannot be ruled out. 



There is other observational support for an anti-correlation between the Mbh — cr* rela- 
tion and Eddington ratio from observations of narrow-line Seyfert 1 galaxies (NLSls). These 
active galaxies typically have relatively small measured BH masses and high accretion rates, 
which results in large Eddington ratios . Seve ral studies (e.g. IWandelll2002l ; iGrupe &: Mathur 
20041 ; iBian fc Zhaoll2004l : iNelson et all l2004h have shown that NLSls tend to lie below the 
Mbh — relation for inactive galaxies. If NLSls are an early stage of AGN evolution, then 
they will presumably move up the Mbh - 0"* plane over time (BH masses growing), becoming 
broad-line AGNs and eventually inactive galaxies harboring supermassive black holes. If the 
objects evolve in this way toward the z = Mbh — relation, the black holes would have 
to be growing more rapidly than the bulges are growing. To explain the tight relationship 
between the final BH and spheroid masses, the gas mass that builds the black hole may 
somehow be set at a given formation stage of the stellar spheroid. It would then be the fate 
of the black hole to grow to a particular mass, determined very early in its developmen t 
by the mass of the spheroid (e.g. iMiralda-Escude fc Kollmeierll2005l ; iMerritt fc Poonll2004l ). 
Feed back from the central engine surrounding the BH may regulate the BH growth pr ocess 



^e.g. iRobertson et al.ll2006l : iBegelman fc Nathll2005l : iDi Matteo et al.ll2005l : Kind 120031 ). but 
this feedback may have little to do with the development of the spheroid. It is also pos- 
sible that NLSls are not a representative class of AGN, and their evolutionary tracks in 
the Mbh — < y * plane are different than those of other types. Observations of AGNs with 
luminosities comparable to NLSls, but at higher redshifts than can be obtained with the 
SDSS, would help clarify the issue. 

Differences in the Mbh — o"* relation may also be due to disk orientation effects, which 
could produce a wide range of incorrect estimates for black hole masses. Given that all BH 
masses were determined by using the same value of the / factor (equation [TJ, the objects 
with high L ho] / L FM may be those in which BH ma sses are underestimated because their true 
/ values are larger. iMclure fc Dunlopl (120011 . 120021 ) argue that a flattened-disk-like geometry 
in the BLR is favored over randomly-oriented orbits, and a flattened geometry viewed over a 
range of orientat ion angles can easil y result in virial BH masses underestimated by a factor of 
three. However, I Collin et al.l (120061 ) showed that the inclination effects are actually minimal 
when using the line dispersions as opposed to the FWHM of the line profile. Although the 
factor of 3 difference in the BH masses (0.5 dex) can explain the apparent dependence of 
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the Mbh — o"* relation on Lboi/^Edd, we cannot settle this issue here because of our poor 
understanding of the BLR geometry. 

Our BH masses were estimated by equati on (|2J), where FWHM was used to characterize 
the line width (IVestergaard fc Peterson 1120061 )). However, some studies showed that the line 
dispe rsion is a less biased parameter in gener al than FWHM for black hole mass estimation 



e.g. 



Peterson et al.l 12004 ; ICoUin et al.ll2006l ). The line di spersion will also minimize the 



effect that may be influenced by relative source inclination (jCollin et al.ll2006l ). 



The SDSS spectroscopic data set has provided a large sample of homogeneously selected 
AGN from which both nuclear and host properties can be measured. The moderate resolution 
and S/N of the SDSS spectra limit the ranges of accessible luminosities, black hole masses, 
and Eddington ratios. However, in principle, the eigenspectrum decomposition technique 
can be applied to any spectrum with adequately strong nuclear and host components, so it 
is possible with other data sets to extend the results of this study to a larger dynamic range. 
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Fig. 1. — Sample object apparent i band magnitudes mi vs. the spectroscopic S/N averaged 
over the % band. The median S/N is 21.6. 
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Fig. 2. — Sample object redshift vs. the spectroscopic S/N averaged over the i band. 




Fig. 3. — AGN monochromatic continuum luminosity at a rest wavelength of 5100A, L 510 o, as 
a function of redshift for the 960 objects in our sample. The AGN luminosity was measured 
after host galaxy subtraction; the measurement of this quantity is discussed in § 13.11 
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Fig. A.— The SDSS spectrum of SDSS J104252.94+041441.1, with z = 0.053, a typical 
example of a high-S/N object in our sample, with S/N = 37.3 (the spectra have a resolution 
of AA/A ~1850). The middle components are the reconstructed AGN and host galaxy 
spectra, the top histogram is the original active galaxy spectrum, and the bottom histogram 
is the fit residual. The gray curve is the modeled active galaxy spectrum, which is red in the 
electronic version. 
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Fig. 5.— The SDSS spectrum of SDSS J212348.61+105348.2, with z = 0.087, a typical 
example of an object with average S/N in our sample, with S/N = 21.7. The middle 
components are the reconstructed AGN and host galaxy spectra with an offset for clarity, 
the top histogram is the original active galaxy spectrum with an offset for clarity, and the 
bottom histogram is the fit residual. The gray curve is the modeled active galaxy spectrum, 
which is red in the electronic version. 
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Fig. 6.— The SDSS spectrum of SDSS J094813.80+401325.9, with z = 0.128, a typical 
example of a low-S/N object in our sample, with S/N = 15.1. The middle components are 
the reconstructed AGN and host galaxy spectra with an offset for clarity, the top histogram 
is the original active galaxy spectrum with an offset for clarity, and the bottom histogram is 
the fit residual. The gray curve is the modeled active galaxy spectrum, which is red in the 
electronic version. 
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Fig. 7. — Example continuum and Fell template fit, for the object described in Figure HI 
The top spectrum is the host-subtracted AGN spectrum with an offset for clarity, the middle 
lines are the featureless continuum and the iron template superposed on the continuum with 
an offset for clarity, and the bottom histogram is the AGN spectrum after subtracting the 
continuum and iron template. The four short horizontal solid lines denote the end points 
used for continuum fitting. The two horizontal dashed lines show the two continuum regions 
used to fit the broken power law. In the electronic version, the green curve is the modeled 
continuum and the red curve is the iron template superposed on the continuum. 



-40- 




4000 4500 5000 5500 6000 6500 7000 

Rest Wavelength (A) 



Fig. 8. — Example continuum and Fell template fit, for the object described in Figure [5j 
The top spectrum is the host-subtracted AGN spectrum with an offset for clarity, the middle 
lines are the featureless continuum and the iron template superposed on the continuum with 
an offset for clarity, and the bottom histogram is the AGN spectrum after subtracting the 
continuum and iron template. The four short horizontal solid lines denote the end points 
used for continuum fitting. The two horizontal dashed lines show the continuum regions 
used to fit the broken power law. In the electronic version, the green curve is the modeled 
continuum and the red curve is the iron template superposed on the continuum. 



-41 - 



50 E 




4000 4500 5000 5500 6000 6500 7000 

Rest Wavelength (A) 



Fig. 9. — Example continuum and Fell template fit, for the object described in Figure EJ 
The top spectrum is the host-subtracted AGN spectrum with an offset for clarity, the middle 
lines are the featureless continuum and the iron template superposed on the continuum with 
an offset for clarity, and the bottom histogram is the AGN spectrum after subtracting the 
continuum and iron template. The four short horizontal solid lines denote the end points 
used for continuum fitting. The two horizontal dashed lines show the continuum regions 
used to fit the broken power law. In the electronic version, the green curve is the modeled 
continuum and the red curve is the iron template superposed on the continuum. 
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Fig. 10.— Distribution of the O III doublet flux ratio (F(5007A)/F(4959A)). The distribu- 
tion shows that most O III doublet ratios are close to the theoretical value of 3:1. 
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Fig. 11. — Example line fit for the object described in Figs. H] and [71 Left panel: Fits to the 
H/3+0 HIAA 4959, 5007A region. The data are shown by the histogram, the fitted broad and 
narrow components are shown as dotted lines, and the sum of all components is shown as a 
solid line. The bottom histogram is the fit residual with an offset of -10 for clarity. Right 
panel: Same as the figure to the left, but for the Ha+[NII]AA 6548,6583A complex. 
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Fig. 12. — Correlation between Ha and H/5 line widths. The solid line gives the BCES 
bisector fit to 59 7 objects. The dash ed line denotes FWHM ffa =FWHM ff)3 . The filled dots 
are the data from lKaspi et al.l (120001 ) . which are the mean FWHM values. A typical FWHM 
error bar for an individual measurement, determined by the mean of all the formal errors of 
FWHM measurements, is shown in the top left part of the panel. 
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Fig. 13. — Distribution of AGN monochromatic luminosity at 5100A, L 5100 . The solid line 
shows the luminosity distribution after subtraction of the host galaxy component, and the 
dotted line shows the distribution without subtraction of the host galaxy. 
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Fig. 14. — The distributions of the Eddington ratio, log(Lboi/-^Edd) (solid histogram, units 
labeled at bottom of plot) and the BH mass, M-q-^/Mq (dotted histogram, units labeled at top 
of plot). The AGN bolometric luminosity, L ho \, was estimated according to the approximate 
relation L bo i « 10AL A (5100A). 
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Fig. 15. — The distribution of the BH masses with the estimated Eddington ratio Lboi/^Edd- 
The apparent anti-correlation is due to parameter interdependencies, and sample selection 
effects (see § 13.51) . The right top horizontal line shows the typical uncertainty of the estimated 
Eddington ratio Lhni / -^rhh , due to the 25% standard deviation in the bolometric correction 
(IRichards et al.ll2006l ) and the BH mass estimation error. The right top vertical line shows 
the statistical BH masses uncertainties. 
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Fig. 16. — An example host galaxy fit to measure the velocity dispersion. The top line is 
the decomposed host galaxy spectrum with an offset 10 for clarity, the middle solid line is 
the fit, and the bottom solid line is the fit residual. The boxes show the masked regions. 
The two large boxes are masks for H (3 and H a, the left and right narrow boxes are masks 
for H7 and [Ol], and the middle box is a mask for bad pixels. 
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Fig. 17. — The relative systematic error and random error in the measurement of velocity 
dispersion, found from simulations. The dotted line shows the relative error between the 
true velocity and the measured noiseless velocity. The dashed line shows the relative error 
between the true velocity and the mean velocity in 100 simulations with noise added, for the 
case of spectroscopic S/N = 15, and host galaxy flux fraction F H = 50%. The solid curves 
show the standard deviations of the relative velocity measurements in 100 simulations. 
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Fig. 18. — The measurement error of the velocity dispersion in the host galaxy sample, 
returned from the vdispfit routine, as a function of host galaxy spectroscopic S/N. 
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Fig. 19. — Same as Fig.[TT], but for the worst case of spectroscopic S/N = 15, and host 
galaxy flux fraction Fjj = 20%. 
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Fig. 20. — Relation between host galaxy g band absolute magnitudes M g , and stellar velocity- 
dispersion o*. Filled circles show the mean values of log a* in sets of objects that span small 
ranges in absolute magnitude; each set contains 50 host galaxies. Vertical error bars show 
the M g bin size, and horizontal error bars show the rms scatter around the mean values of 
a*. The solid line shows the best fit linear relation (in logarithmic space). 
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Fig. 21. — Host galaxy absolute g band magnitude vs. BH mass. The data are binned by 
M g in the same way as in Fig.[20j Filled circles show the mean value of M-q^/Mq in each 
bin. The horizontal error bars show the rms scatter around the mean BH values, and the 
vertical error bars show the M g bin size. 
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Fig. 22. — Black hole mass vs. stellar velocity dis persion for individua l points. The dot- 
dashed line shows the best fit for inactive galaxies (jTremaine et al.ll2002l ) 
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Fig. 23. — Black hole mass vs. stellar velocity dispersion in binned data. The solid line is 
the best-fitting linear relation (in logarithmic space), which gi v es the slope 3.34±0.24 and 
intercept 7.92 ±0.02. The dashed line is the iTremaine et al.l (120021 ) relation for inactive 
galaxies. Each point represents 50 objects which have been binned by host galaxy absolute 
magnitude, as in Figs. [20] and ED The crosses show the mean value in each bin, and error 
bars are the standard deviation of the mean in each bin. 
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Fig. 24. — Black hole mass vs. stellar velocity dispersion for several different data sets, 
covering a larger dynamic range t han Fig.l23l The fille d circ les are the da t a from this paper. 
The asterisks are the data from lOnken et al.l (120041 ) and iNelson et al.l (120041) . who used 



reverbe ration mapping measurements. The open s quares are the data f rom I Greene fc Ho 



(2006a), and the open triangles are the data from iPeterson et al.l (120051 ) and iBarth et al. 
J2004h for NGC 4395 and POX 52 respectively. The encircled objects are N LSls. All the 
l iteratu re data values and their uncertainties are available in the table from I Greene fc Ho 
(j2006al ). The solid line is the best-fit relation for all of the data (this paper and the literature 
data combined). For comparison, the dash ed line shows the best-fit relation for inactive 
galaxies, as found by lTremaine et al.l (120021 ). Typical uncertainties in the measurements of 
BH mass and velocity dispersion for individual objects are shown in the bottom right corner. 
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Fig. 25. — The Mbh — c* relationship for two samples divided by the mean value of the 
Eddington ratio, L bo i/L Edd . Mean values of M B h and a* for objects with L bo i/L Edd below 
the median value of 0.027, are shown with filled squares; open triangles show the mean values 
for objects with Lboi/-^Edd above the median value. Each sample has contributions from 451 
objects, and each point represents the mean values for at least 30 objects. The error bars are 
the standard deviation of the mean i n each bin. The s olid li ne shows the best fit relation from 
equation (jSJ . The dashed line is the iTremaine et al.l (120021 ) relation for inactive galaxies. 
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Fig. 26. — The BH mass difference between the measured BH mass and BH mass predicted 
by equation (JBJ) versus the Eddington ratio, Lboi/^Edd- 
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Fig. 27. — The distribution of the host galaxy velocity dispersion with the AGN monochro- 
matic luminosity at 5100A. 
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Fig. 28. — The distribution of the host galaxy g band luminosity L H , with the AGN 
monochromatic luminosity at 5100A. The apparent correlation may be due in part to sample 
selection effects. 
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Fig. 29. — The Mbh — cr* relationship in different redshift ranges. Objects with redshifts 
from to 0.1 are represented with open triangles, which show the weighted average values 
in bins containing measurements from at least 25 objects. Filled squares show the mean 
values for objects with redshifts from 0.1 to 0.2; each point has contributions from at least 
40 objects. Crosses show the mean values for objects with redshifts from 0.2 to 0.452; each 
point represents data from at least 30 objects. The error bars are the standard deviation of 
the mean in each bin. The solid line shows the best fit relation from equation (jSj). With the 
slope fixed at the value in equation®, the dotted, dashed and dot-dashed lines show the 
best fits for the three bins with redshifts from low to high respectively. 



Table 1. Sample Properties 



Name 


z 


mi F H 4> 


S/N 


S/AT(Host) 


Mg(HOSt) 1 


og I/ 5 ioo (1) 1 


og L 5 ioo(2) 






FWHM H/3 


FWHM Ha 


log M B h 


ibol/^Edd 


(1) 


(2) 


(3) (4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


000611.54+145357.2 0.1186 15.79 0.41 2.30 


25.9 


15.7 


-22.47 


39.52 


39.90 


181.1+13.4 


188.9 


3.39+0.11 


3.13+0.05 


6.96+ 0.07 


0.1568 


000729.98-005428.0 


0.1454 


17.55 0.53 3.54 


15.3 


9.6 


-20.71 


39.08 


39.65 


121.2+17.8 


127.5 


3.30+0.32 


2.92+0.07 


6.63+ 0.19 


0.1217 


000805.62+145023.3 


0.0455 


14.93 0.65 3.70 


44.9 


29.2 


-20.78 


38.67 


39.37 


173.1+ 7.9 


173.4 


6.09+0.25 


7.02+0.07 


6.88+ 0.09 


0.0266 


000813.22-005753.3 


0.1393 


16.45 0.46 3.64 


27.5 


13.3 


-21.75 


39.61 


39.97 


221.8+17.6 232.9 


5.21+0.24 


3.24+0.05 


7.39+ 0.10 


0.0717 



to 

Note. — Col.(l):Name (SDSS J). Col. (2): Redshift. Col. (3): i band magnitude. Col. (4): Host galaxy fraction. Col. (5): Galaxy classification angle (degrees). Col. (4): 
i band spectroscopic S/N. Col. (7): Host galaxy S/N. Col. (8): Host galaxy absolute g band magnitude. Col. (9): Z/5100 (ergs s _1 A _1 , host subtracted). Col. (10): Z/5100 
(ergs s _1 A — 1 , host not subtracted). Col. (11): Host galaxy velocity dispersion (km s~ 1 ), uncorrected for finite fiber diameter. Col. (12): Host galaxy velocity dispersion 
(km s _1 ), corrected for finite fiber diameter by equation gj. Col. (13): FWHM H(3 (10 3 km s" 1 ). Col.(14): FWHM Ha (10 3 km s" 1 ). See §|jl2]for the error estimates of 
FWHMna and FWHM H/3 . Col. (15): M B h (Mq), calculated from equation (0 by FWHM H(3 . The symbol '*' indicated where FWHM H(3 was obtained from FWHM ffo 
by equation {3). The A/bh uncertainties are formal uncertainties, actual uncertainties are probably dominated by systematics in BLR geometry. Col. (16): ibol/^Edd- 
From Col. (13) to Col. (16), some quantities are set to zeros. These default values are given for entries with bad spectra or measurements (see § 13.31 . Table 1 is available 
in the electronic edition. A portion is shown here for guidance regarding its form and content. 
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Table 2. Mean properties in host galaxy luminosity bins 



Mg(H) 


Mgj a i n t(il) 




logM BH 


A log Mbh 


log cr* 


A log (7* 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




99 89 


93 RR 


R 31 


n ^1 

U. O -L 


9 38 
z. oo 


n in 

U . 1 u 


- ZZ .UU 


99 ^ 
-zz. oo 


99 81 

- ZZ. O 1 


RIO 

. 1 u 


n 3fi 


9 3fi 
z. ou 


n 1 1 

U . J. -L 


99 41 


99 31 

"ZZ.Ol 


99 ^ 


7 Q9 
i .yz 


n 3& 

u.oo 


9 33 
z.oo 


n 1 3 

u. io 


99 99 
-zz .zz 


99 14 


99 31 


7 Qfi 


u.oo 


9 30 
z. ou 


nil 

U. X X 


99 OR 
- ZZ .uu 


99 nn 


99 1 3 
-zz. ±o 


7 Q1 


u.ou 


9 39 
z. oz 


n 1 1 

U . 1 X 


91 Q 1 ^ 
-z i .yo 


91 Q1 

-Z l.jl 


91 QQ 

-Z 1 . c/c/ 


7 QD 




9 97 
z.z / 


19 

U . 1 z 


91 &zl 

_ Z 1 .O'i 


91 7 ^ 

-Z ±. ( u 


91 Q1 
-z i . y i 


7 7zL 


zL9 


9 9^ 
z. zo 


PI DQ 

U.Uc/ 


-21.70 


-21.64 


-21.76 


7.57 


0.38 


2.22 


0.10 


-21.59 


-21.54 


-21.64 


7.63 


0.34 


2.25 


0.12 


-21.49 


-21.45 


-21.54 


7.62 


0.40 


2.24 


0.09 


-21.37 


-21.30 


-21.44 


7.60 


0.48 


2.23 


0.12 


-21.23 


-21.16 


-21.30 


7.70 


0.33 


2.18 


0.10 


-21.08 


-21.02 


-21.16 


7.53 


0.28 


2.20 


0.11 


-20.95 


-20.89 


-21.01 


7.44 


0.31 


2.17 


0.12 


-20.79 


-20.71 


-20.89 


7.50 


0.35 


2.16 


0.12 


-20.59 


-20.46 


-20.70 


7.38 


0.45 


2.15 


0.10 


-20.29 


-20.08 


-20.46 


7.48 


0.33 


2.11 


0.11 


-19.60 


-18.41 


-20.07 


7.20 


0.35 


2.11 


0.14 



Note. - 903 objects with reliable parameter measurements were dividend into 
18 bins; each bin contains 50 objects except the last bin which contains 53 ob- 
jects. Col.(l): Mean host galaxy g band absolute magnitude in the bin. Col. (2): 
Faintest host galaxy g band absolute magnitude in the bin. Col. (3): Brightest 
host galaxy g band absolute magnitude in the bin. Col. (4): The weighted mean of 
the log Mbh (M ) in the bin. Col. (5): The dispersion of log Mbh (M ) in the bin. 
Col. (6): The weighted mean of velocity dispersion logo - * in the bin. Col. (7): The 
dispersion of logo - * in the bin. 
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Table 3. Partial correlation analysis for the dependence of AM BH on Eddington ratio 



X 


r L bo i/L Edd ,x 


r AM,x 


fAM,L bol /L Edd -x 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


-^5100 


0.377 


-0.081 


-0.591 


< 0.01% 


FWHM 


-0.882 


0.577 


-0.176 


< 0.01% 




0.023 


-0.563 


-0.682 


< 0.01% 




0.030 


-0.044 


-0.576 


< 0.01% 


z 


0.137 


0.005 


-0.583 


< 0.01% 



Note. - Partial correlation coefficients for AMbh with 
Eddington ratio. Col.(l): The control parameters in the 
partial correlation test. Col. (2): The correlation between 
Eddington ratio and the control parameter. Col. (3): The 
correlation between AMbh and the control parameter. 
Col. (4): The partial correlation coefficient for the inverse 
correlation of AM BH with Eddington ratio when the influ- 
ence of the control parameter is accounted for. Col. (5): The 
significance of the partial correlation. 
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Table 4. Partial correlation analysis for the dependence of L5100 on °"* 



X 








P 


(1) 


(2) 


(3) 


(4) 


(5) 


L H 


0.634 


0.533 


0.354 


< 0.01% 


z 


0.802 


0.524 


0.292 


< 0.01% 


All a 






0.262 


< 0.01% 



Note. - Partial correlation coefficients for 
AM BH with Eddington ratio. Col.(l): The 
control parameters in the partial correlation 
test. Col. (2): The correlation between L 510 o 
and the control parameter. Col. (3): The cor- 
relation between a* and the control parameter. 
Col. (4): The partial correlation coefficient for 
the control parameter. Col. (5): The signifi- 
cance of the partial correlation. 

a Second-order partial correlation test, ac- 
counting for both control parameters. 
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Table 5. Partial correlation analysis for the dependence of AMbh on redshift 



X 


T z,x 


T'AM,x 


r AM,z;x 


P 


(1) 


(2) 


(3) 


(4) 


(5) 


-^5100 


0.802 


-0.081 


0.118 


< 0.04% 


FWHM 


0.257 


0.577 


-0.181 


< 0.01% 




0.524 


-0.563 


0.428 


< 0.01% 


L H 


0.707 


-0.044 


0.052 


11.8% 


Lbol/L>Edd 


0.142 


-0.576 


0.108 


0.11% 


All a 






0.059 


7.64% 



Note. - Partial correlation coefficients for 
AM B h with redshift. Col.(l): The control pa- 
rameters in the partial correlation test. Col. (2): 
The correlation between redshift and the con- 
trol parameter. Col. (3): The correlation between 
AMbh and the control parameter. Col. (4): The 
partial correlation coefficient for the control pa- 
rameter. Col. (5): The significance of the partial 
correlation. 

a Fifth-order partial correlation test, accounting 
for all 5 control parameters . 
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